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Abstract 

The correlated behavior of electrons determines the structure and optical properties of molecu- 
les, semiconductor and other systems. Valuable information on these correlations is provided by 
measuring the response to femtosecond laser pulses, which probe the very short time period during 
which the excited particles remain correlated. The interpretation of four-wave-mixing techniques, 
commonly used to study the energy levels and dynamics of many-electron systems, is complicated by 
many competing effects and overlapping resonances. Here we propose a coherent optical technique, 
specifically designed to provide a background-free probe for electronic correlations in many-electron 
systems. The proposed signal pulse is generated only when the electrons are correlated, which 
gives rise to an extraordinary sensitivity. The peak pattern in two-dimensional plots, obtained 
by displaying the signal vs. two frequencies conjugated to two pulse delays, provides a direct 
visualization and specific signatures of the many-electron wavef unctions. 

PACS numbers: 71.35.Cc, 73.21.La, 78.47.+p 
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Predicting the energies and wavefunctions of interacting electrons lies at the heart of 
our understanding of all structural, optical, and transport properties of molecules and 
materialgii^i'^'^i^i^i^i^i'i^. The Hartree-Fock (HF) approximation provides the simplest de- 
scription of interacting fermions^'^. At this level of theory each electron moves in the aver- 
age field created by the others. This provides a numerically tractable, uncorrelated-particle 
picture for the electrons, that approximates many systems well and provides a convenient 
basis for higher- level descriptions. Electronic dynamics is described in terms of orbitals, 
one electron at a time. Correlated n-electron wavefunctions, in contrast, live in a high (3n) 
dimensional space and may not be readily visualized. Deviations from the uncorrelated 
picture (correlations), are responsible for many important effects. Correlation energies are 
comparable in magnitude to chemical bonding energies and are thus crucial for predicting 
molecular geometries and reaction barriers and rates with chemical accuracy. These energies 
can be computed for molecules by employing a broad arsenal of computational techniques 
such as perturbative corrections^^, configuration interaction^, multideterminant techniques^, 
coupled cluster theory^ and time dependent density functional theory TDDFT-"^*^. Cor- 
relation effects are essential in superconductors^i^^ and can be manipulated in artificial 
semiconductor nanostructure9^>^ii^. The fields of quantum computing and information are 
based on manipulating correlations between spatially separated systems, this is known as 
ent anglemenlr^ . 

In this article we propose a nonlinear optical signal that provides a unique probe for 
electron correlations. The technique uses a sequence of three optical pulses with wavevec- 
tors ki,k2 and ks, and detects the four wave mixing signal generated in the direction 
k5= ki+k2— ks by mixing it with a fourth pulse (heterodyne detection) .-i^ We show that this 
correlation- induced signal Sci(t3,t2,ti), which depends parametrically on the consequentive 
delays ti, t2, ts between pulses, vanishes for uncorrelated systems, providing a unique indi- 
cator of electron correlations. This technique opens up new avenues for probing correlation 
effects by coherent ultrafast spectroscopy. 

Starting with the HF ground-state (g) of the system, each interaction with the laser 
fields can only move a single electron from an occupied to an unoccupied orbital. The 
first interaction generates a manifold (e) of single electron-hole (e-h) pair states. A second 
interaction can either bring the system back to the ground state or create a second e-h pair. 
We shall denote the manifold of doubly excited states as / (Fig. 1). We can go on to generate 
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manifolds of higher levels. However, this will not be necessary for the present technique. 
The quantum pathways (i) and (ii) contributing to this signal can be represented by the 
Feynman diagrams^^ shown in Fig. 1. Each diagram shows the sequence of interactions of 
the system with the various fields and the state of the electron density matrix during each 
delay period. We shall display the signal as Sci{^3,^2,ti), where ^3, Q2 are frequency 
variables conjugate to the delays ^3 and t2 (Fig. 1) by a Fourier transform 

/•oo /»oo 

Sci{^3,^2,tl) = / dt2dtsSci{t3,t2,ti)exp{iQ2t2 + i^3t3), 
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with ti fixed. This yields an expression for the exact response function 

Sci{n3,n2,h = o) (1) 

f^gef^ef fe' f^e' g f^ge' f^e' f f^fef^eg 

where for simplicity we set ti = 0. Two-dimensional correlation plots of fl2 vs. fl3 then 
reveal a characteristic peak pattern, which spans the spectral region permitted by the pulse 
bandwidths. The two terms in the brackets correspond respectively to diagrams (i) and 
(ii) of Fig. 1. Here fi^^' are the transition dipoles and uji,,^/ are the transition energies 
between electronic states, shifted by the pulse carrier frequency Uq, i.e., u^'g = e^' — eg —ujq, 
ojfe = ej — ee — ujQ and Ufg = ef — eg — 2u!o. This shift eliminates the high optical frequencies. 
The carrier frequency of the three beams, uq, is held fixed and used to select the desired 
spectral region. In Eq. ([1]) we have invoked the rotating wave approximation (RWA), and 
only retained the dominant terms where all fields are resonant with an electronic transition. 

In both diagrams, during t2 the system is in a coherent superposition (coherence) between 
the doubly excited state / and the ground state g. This gives the common prefactor (^2 — 
ujfg)~^ . As is scanned, the signal will thus show resonances corresponding to the different 
doubly-excited states /. However, the projection along the other axis (^3) is different for 
the two diagrams. In diagram (i) the system is in a coherence between e' and g during 
^3. As ^3 is scanned, the first term in brackets reveals single excitation resonances when 
0.3 = uJe'g. For the second diagram (ii) the system is in a coherence between / and e' during 
^3. This gives resonances at = ujje in the second term in the brackets. Many new peaks 
corresponding to all possible transitions between doubly and singly excited states ojfe should 
then show up. 
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The remarkable point that makes this technique so powerful is that the two terms in 
Eq. ([1]) interfere in a very special way. For independent electrons, where correlations are 
totally absent, the two e-h pair state / is simply given by a direct product of the single 
pair states e and e', and the double-excitation energy is the sum of the single-excitation 
energies = eg + eg', so that uo^'g = ^fe and the two terms in the brackets exactly cancel. 
Density functional theoryii^iii, when implemented using the Kohn-Sham approach, gives a 
set of orbitals that carry some information about correlations in the ground state. The 
signal calculated using transitions between Kohn-Sham orbitals will vanish as well. This 
will be the case for any uncorrelated-particle calculation that uses transitions between fixed 
orbitals, no matter how sophisticated was the procedure used to compute these orbitals. We 
expect the resonance pattern of the two-dimensional Sci signal to provide a characteristic 
fingerprint for electron correlations. 

The following simulations carried out for simple model systems which contain a few 
orbitals and electrons, illustrate the power of the proposed technique. Doubly excited 
states can be expressed as superpositions of products of two e-h pair states. Along ^3 we 
should see the various doubly excited states at uofg^ whereas along we observe the various 
projections of the / state onto single pair states uj^'g and the differences ojje = ojfg—ojeg- The 
2D spectra thus provide direct information about the nature of the many body wavefunctions 
that is very difficult to measure by any other means. The patterns predicted by different 
levels of electronic structure simulations provide a direct means for comparing their accuracy. 
We used a tight-binding Hamiltonian H = Hq + He + H^. The single-particle contribution 
Hq contains orbital energies and hoppings 



where c„j and dn^ are electron and hole annihilation operators respectively, and the sum- 
mations run over spin-orbitals. We assume equal hopping t for electrons and holes. The 
many-body term responsible for correlations 



contains only direct Coulomb couplings. The electron-electron, hole-hole and electron-hole 
interaction are denoted V^^, yhh^ V^^, respectively. Values of t and Coulomb integrals 
Vqq, Vqi (subscripts and 1 indicate the sites) were derived by fitting emission spectra of 
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coupled quantum dots^^. Owing to the nature of quantum dots states, we can assume that 
yeh yee ^ yhh is -yy^ pj^oose V^'^ = V^^ = 1.2V^^' and usc thesc values for all orbitals. 
Hi describes the dipole interaction with the laser pulses, = —E(t)fijnim2d"m,2Cmi + h.c, 
where E{t) is the light field and fimim2 ^i-re local dipole moments of various orbitals mi, m2- 
Uq was tuned to the single-particle optical gap energy (1 3 eV). 

Even though our parameters are fitted to QDs, the overall picture emerging from the 
calculations can be applied to a wider class of systems, whose optical response is determined 
by correlated e-h pairs. We have employed an equation of motion approach for computing 
the signal. Many-body states are never calculated explicitly in this algorithm. Instead, we 
obtain the signal directly by solving the Nonlinear Exciton Equations (NEE)^i2i. These 
equations describe the coupled dynamics of two types of variables representing single e-h 
pairs: Bm = {dm2Cmi) (here m = {1711,1712) stands for both the electron index mi and the 
hole index 1112) and two pairs Ymn = {dm2Cmidn2Cni) ■ For our model the NEE is equivalent 
to full CI, and yield the exact signal with all correlation effects fully included. This signal 
provides a direct experimental test for many-body theories, which use various degrees of 
approximations to treat electron correlations. We compare the exact calculation (NEE) 
with the time dependent Hartree Fock (TDHF) theory, which is an approximate, widely used 
technique for treating correlations by factorizing the Y variables into {dm2Cmi) {dn2Cni) — 
(dnjCrni) {dm2Cni)- This assumcs that two e-h pairs are independent and we only need to 
solve the equations for ((i„2C„^). Correlation within e-h pairs is nevertheless retained by this 
level of theory, as evidenced by the finite Sci signal. The equations of motion derived using 
both levels of factorization are solved analytically, yielding the exact and the TDHF Sci 
signals. The TDHF solutions have the following structure: a set of single-particle excitations 
with energies and the corresponding transition dipole moments are obtained by solving 
the linearized TDHF equations. Many-particle state energies are given by sums of these 
elementary energies. Two-particle energies are of the form ea + 6/3. This approximation is 
the price we pay for the enormous simplicity and convenience of TDHF. Correlated many- 
electron energies computed by higher level techniques do not possess this additivity property. 

In general the TDHF signal contains a different number of resonances along Q2 than the 
exact one. Their positions, Ufg = ef — eg, also differ since the former uses the additive 
approximation for the energies ef. The ^3 value of each resonance in the exact simulation 
is given by either cjg'g (first term in brackets in Eq. [T]) or Ufg (second term in the brackets). 
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The simulations of the TDHF response function presented below show fewer peaks than in 
the exact calculation. This dramatic effect reflects direct signatures of the correlated two 
e-h pair wavefunction, which are only revealed by the Sci technique. 

We first consider a simple model, consisting of a single site with one valence orbital and 
one conduction orbital (Fig. 2A). The energy of the (spin-degenerate) single-pair state is 
ee = —V^^. The only two-pair state has energy e/ = —AV^^ + V^'^ + V^^, compared with 
e/ = 2ee in the TDHF approximation (the TDHF double-excited energies and frequencies will 
be marked with a bar: e/, ojfe)- Thus, the exact signal has two peaks at {VL^, VL2) = iyJegi ^fg) 
and (ufeyUJfg), while TDHF predicts only one peak at {Q3,Q2) = iy^eg.^fg) = {^fe,^^fg)- 

This high sensitivity to correlation effects is general and is maintained in more complex 
systems. In Fig. 2B we consider a system with two valence orbitals with a splitting A and 
one conduction orbital. It has 2 single e-h pair transitions ei, 62 with energies ei = —V^^ 
and 62 = —V'^^ + A. The exact spectrum contains 8 peaks, with e/ energies being sums of 
all quasiparticle interactions and hole level energies: e/i = —AV^^ + V^'^ + V^^, ej^ = ej^ + A, 
ef.^ = ef^ +2 A. Within TDHF we find e/^ = 2ei, e/j = ei -|-e2, e/3 = 2e2 and the 2D spectrum 
only shows 4 peaks. 

The simple energy- level structure of the two systems (Fig. 2A,B), whereby COf^ = 
ujegi allows an insight into the differences in predictions of the two response func- 
tions. We recast Eq. ([1]) for the exact signal in a slightly different form: 
[0.2 — ^fg)~~^ (^3 — ^e'gY^ (^3 — ^/e) ""^ (^H dipolcs for this system are equal /igg' = Z^ge = 
/ie'/ etc.). The corresponding TDHF expression is: (^2 — ^fg) ^ (^3 — ^e'gY'^ ■ The different 

dependencies reflect different numbers of resonances with different FWHM along the VL2, 
axis. The double resonance in TDHF is split into two resonances in the exact expression. 

Fig. 2C shows the signal from two coupled quantum dots, each hosting one valence and 
one conduction orbital. The Sci signal contains a rich peak structure, reflecting the 4 
(10) many body levels in the single- (double-) excited manifold. Again, the TDHF method 
misses many peaks. In this case, unlike the two previous systems, TDHF does not show all 
possible resonances along the 1^2 axis. This is because one of the single-excited states (ei) 
is not optically allowed. In TDHF any / state, constructed as a direct product of ei with 
another state Cj {i = 1, ... ,4) is forbidden, thus we have only 6 resonances. In the exact 
calculation, the / states are not direct products, so we see all possible resonances along fl2- 
The differences between the TDHF and exact spectra in all these examples, illustrate the 
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sensitivity of the proposed signal to the correlated wavefunction. 

Computing electron correlation effects, which are neglected by HF theory, consti- 
tutes a formidable challenge of many-body theory. Each higher-level theory for electron 
correlationsii^ is expected to predict a distinct two-dimensional signal, which will reflect 
the accuracy of its energies and many-body wavefunctions. The proposed technique thus 
offers a direct experimental test for the accuracy of the energies as well as the many-body 
wavefunctions calculated by different approaches. Time dependent density functional theory 
TDDFT within the adiabatic approximation extends TDHF to better include exchange and 
correlation effects^*^. However, the two are formally equivalent and yield a similar excited- 
state structure^^. The two-dimensional peak pattern of TDDFT will suffer from the same 
limitations of TDHF. 

We can summarize our findings as follows: at the HF level which assumes independent 
electrons, the Sci signal vanishes due to interference. TDHF (or TDDFT) goes one step 
further and provides a picture of independent transitions (quasiparticles). Here the signal 
no longer vanishes, but shows a limited number of peaks. When correlation effects are 
fully incorporated, the many-electron wavefunctions become superpositions of states with 
different numbers and types of e-h pairs. The Q2 and ^3 axes will then contain many more 
peaks corresponding to all many body states (in the frequency range spanned by the pulse 
bandwidths), which project into the doubly-excited states. Thus, along ^2 the peaks will be 
shifted, reflecting the level of theory used to describe electron correlations. Along ^3, the 
effect is even more dramatic and new peaks will show up corresponding to splittings between 
various levels. This highly-resolved two dimensional spectrum provides a invaluable direct 
dynamical probe of electron correlations (both energies and wavefunctions). 

Signals obtained from a similar pulse sequence, calculated for electronic transitions in 
molecular aggregates'^ and molecular vibrations'^ show the role of coupling between Frenkel 
excitons. A conceptually related NMR technique known as double quantum coherence re- 
veals correlation effects among spins. The technique showed unusual sensitivity for weak 
couplings between spatially remote spins and has been used to develop new MRI imaging 
techniques2^i2^. Here we have extended this idea to all many-electron systems. The proposed 
technique should apply to molecules, atoms, quantum dots and highly correlated systems 
such as superconductors. It has been recently demonstrated that two-exciton couplings can 
be controlled in onion-like semiconductor nanoparticles with a core and an outer shell made 
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of different materialsl^,. Nonlinear spectroscopy of tlie kind proposed here could provide 
invaluable insights into the nature of such two-exciton states. 
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Figure Captions 



Fig. 1 Left: many-body states connected by transitions dipoles, including the ground 
state g, the manifold of single e-h pairs e and the manifold of two-pair states /. Right: the 
two Feynman diagrams contributing to the correlation-induced signal Sci (Eq[T]). ti are the 
time delays between laser pulses. For independent electrons ujfe = oJeig and diagrams of type 
(i) and (ii) cancel in pairs. 

Fig. 2 Absolute value of the exact and the TDHF Sci signals for three model systems. 
Energies on the axes are referenced to the carrier frequency, which excites interband tran- 
sitions. Each system has A^(A^ + l)/2 doubly-excited levels, where is the number of 
single-excited levels. (A) 2 orbital system with A^ = 1 {y^^ = 194.4cm~^), (B) 3 orbital 
system with N = 2 {V^^ = 194.4cm~-'^), (C) 2 coupled systems of type (A), but with dif- 
ferent gaps, = 193.9 cm~^ for one dot and 78.1cm~^ for the other and t = 59.2cm~^, 
A^ = 4. The corresponding orbitals and splittings are given schematically below each panel. 
For systems (A) and (B), the TDHF misses half of the resonances along ^3, while for (C) it 
misses 4 out of 10 resonances along 
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